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Abstract. Values for the vacuum energy of scalar fields under Dirichlet and Neumann bound- 
s/3 . ary conditions on an infinite clylindrical surface are found, and they happen to be of opposite 
signs. In contrast with classical works, a complete zeta function regularization scheme is here 
applied. These fields are regarded as interesting both by themselves and as the key to describ- 
ing the electromagnetic (e.m.) case. With their help, the figure for the e.m. Casimir effect in 
the presence of this surface, found by De Raad and Milton, is now confirmed. 
C£ 1 PACS: 03.70.+k, 12.20.-m, 42.50.Lc 
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Introduction. The Casimir effect |1J is caused by zero-point fluctuations of quantum fields 
in the vacuum when they are modified by the introduction of boundaries or constraints. The 
particular features of the resulting force depend on the nature of the field, boundaries and 
boundary conditions (b.c.) imposed (in some cases, an interpretation in terms of the radiation 
pressure has been made). It arises in different areas of physics, including quantum electro- 
dynamics, theory of hadrons, and cosmology. Further, there has been a series of attempts to 
explain the phenomenon of sonoluminiscence on the basis of the Casimir effect, which, although 
unsuccessful (see || and refs. therein) still gives an idea of the importance attributed to its 
role in modern physics. 

Evaluating the zero-point energy is usually an involved problem, and a great deal of 
techniques have been thought up for this purpose. Some reasons have made of zeta function 
regularization ||-[EJ a relatively popular method. Particularly, the variant we employ in this 
paper has similarities with the technique developed in ref.JTJ for a cosmological problem. In 
ref. 0, it was successfully applied to the cases of the sphere and the circle reobtaining known 
results [HJ-flnj] and finding some new ones. Old as the spherical problem may look, it has 
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been repeatedly revisited in such recent papers as refs. |T5], [TjJ. However, as far as we know, 
a complete zeta function approach has never been applied to the problem of the cylinder, first 
studied in ref . ||L7H . This is precisely what we plan to do in the present work. 

The spectral zeta functions, made from some operator's eigenmode set {uj n } = A4, will 
be called 

( M (s)=^n S - (1) 



Quite often, it is even more convenient to use the dimensionless version (m (s) — ( — ) , 

M n V P J 

where fi is an arbitrary scale with mass dimensions. In most cases of interest, we have to 
consider spectra of the same sort as that of a free particle, and u n grows with n without bound. 
Therefore, strictly speaking, these identities hold only for Re s > Sq, being sq a positive number 
given by the rightmost pole of (m{ s )- Nevertheless, this function admits analytic continuation 
to other values of s, including the negative reals. 

The vacuum energy density e will be obtained by zeta-regularization of the mode sum 
-^c<j„, (we have adopted the typical QFT units, i.e., h — c — 1). If the result found is 

2 n 

not finite yet, one can add a principal part (P) prescription as in ref. 0. Writing everything 
together, 



e{n) 



-/i (m(s); s 

Z M 
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where P means extraction of the principal part. Obviously, for this procedure to be operational, 
we need the analytic continuation of (m( s ) to s = — 1. In our case, as we have an infinite 
cylinder, we consider e as the 'linear' density of energy per length unit. 

The present problem can be related to that of a circle in 2 + 1 dimensions. If {u> n } = C 
denotes the set of all the eigenfrequencies, the complete zeta function for the circle will be 
Cc(z) = 51 UJ n Z - Given that the ones for the cylinder are M. = {Juj% + k 2 , uj n G C,k 6 R}, 
the new zeta function will be 
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where B stands for the Euler beta function. Laurent-expanding around s 
that 

~Cc(-2) 



Cm(s) = - 



s + 1 



+ ln(2) - - Rc(-2) + Cc(-2) + 0(s + 1) 



(3) 

-1 one realizes 
(4) 
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Thus, to determine (m{ s ) at s = — 1 means having to know (c{ z ) at and around z = —2. For 
this reason, it is not possible to directly apply the results in ref . for the circle, which only give 
information about (c(z) at z — —1. The residue of Cm( s ) at s = — 1 is Cc( — 2) itself. In fact, 
it turns out that the sum of the internal and external contributions to this quantity vanishes. 
This can be anticipated by the parity of the Seeley - de Witt coefficient related to the possible 
divergence (actually, the same as the coefficient for the scale anomaly of the one-loop energy) 
and by the fact that we are including the internal and external spaces to an infinitely thin 
boundary in a globally flat manifold of odd space dimension — see e.g. 0, || (pages 107-111), 
fl8|| or refs. therein. Such a vanishing means two things: 

• First, the whole result is finite and the principal part prescription, as well as the use of 
the scale fi, are actually unnecessary. Then, (0) reduces to 

s = \Cm{-1). (5) 

• Second, the finite value which appears at the end depends only on ^(—2), i.e. 

Cm(*) = ^Cc(-2) + 0(3 + 1). (6) 

Z7T 

Therefore, to know only the value of (c at s = — 2 — without knowing the next term in 
its Taylor expansion — would be of no help. 

The present paper is organized as follows: In the next section we consider the Maxwell 
equations and how the e.m. problem decomposes into two situations with just scalar fields. 
Afterwards, the relevant zeta functions are introduced in sec. 3. Then, the energetic contribu- 
tions of the Dirichlet and Neumann modes are separately found in sees, 4 and 5, respectively. 
The total result and related comments appear in sec. 6. 

Solutions of the Maxwell equations for a cylinder. The eigenfrequency set of an elec- 
tromagnetic problem is the one dictated by the Maxwell equations under the b.c. in question. 
Therefore, our first step will be to solve the Maxwell equations. We will look for solutions in 
cylindrical coordinates of the form 

E u Bj = e iuJt e ikz e im8 R{r). (7) 
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We observe that in ([7p u > k, which allows us to boost to a frame in which k = 0. In 
this frame, Maxwells equations read 



im 

B z 

dB z 

dr 
d(rB e 

dr 



= —iuiE r 
= icuEg 
— imB r 



iurE 7 



and 
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— E z 

dE z 
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= iuiB r 
= —iusBq 

— imE r = —iurB z 



(8) 



In fact, (§) are two decoupled systems of three equations each, which can be solved using well 
known techniques. Inside the cylinder, after imposing regularity at r = (which excludes the 
Y m solutions), we find two sets of solutions: The first one is given by 



B z = J m {u>r), 
while the second one reads 

E z = J m {ujr), 
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Or 



B r = Jm(ur) 
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% dJ m (tur) 
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dr 



(9) 



(10) 



where we only list the radial parts of the solutions, and the components not listed vanish. 
Calling n the normal vector to the boundary, the b.c. for an infinite and perfectly conducting 



0, n x E 



0, imply E = E z = 0, B r 



0. 



cylinder of radius a, which read n ■ B 

i — a 

and, thus, J' m {uja) = for @, and J m (ua) = for fllPl) . 

In regard to the solutions outside of the cylinder, we have to bear in mind that the 
Y m functions are no longer to be ruled out. Thus, we may start by replacing J m {ujr) — > 
AJ m {ujr) + BY m {ujr) in @ and (0). Next, the conditions at infinity impose constraints on A 
and B. This question has been addressed in some detail in for the sphere, but the same 
type of calculation may be carried over to the present case. The outcome is that, in practice, 
everything works out as if the actual combination was J m (ujr) + iY m (ujr) = H^(ur) (or its 
conjugate). Once these solutions have been obtained, one just has to boost back to the original 
frame. 

Summing up, the Maxwell equations for cylindrical waves under perfectly conducting 
b.c. on an infinite cylinder of radius a have two sets of solutions, corresponding to the transverse 
electric (TE) and transverse magnetic (TM) modes. TE modes are subject to Dirichlet b.c, 
and TM modes obey Neumann b.c. Thus, when restricting the problem to the circle, the set 
of modes may be written as C = T> U A/", where D stands for the Dirichlet modes and TV for the 
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Neumann modes. These two sets are further subdivided, depending on whether we consider the 
part of the solution inside or outside of the cylinder. In the first case, we add the requirement of 
regularity in the interior, while in the second we have to impose physically suitable conditions 
at infinity. For the internal TE modes (|10D , the u's are such that uoa must be nonvanishing 
zeros of the Bessel function J m , i.e., 



J m {ijja) =0, m G Z, (11) 

where m is the angular momentum number, while the internal TM eigenfrequencies (|9|) obey a 
Neumann condition of the form 

J' m {uja) = 0, m E Z. (12) 

As for the external modes, and in view of the the analytical methods we plan to apply, it is 
formally enough to replace every Bessel function J m by a Hankel function H$ (justified by 
H). In consequence, I m functions should be replaced by K m functions (See below). 

Partial-wave and complete zeta functions. As we have just seen, the e.m. problem can be 
regarded as the sum of two situations involving a massless scalar field: in the first, it is subject 
to flllD , and, in the second, has to satisfy (|12|). Let's now consider an adequate formalism for 
each part. 

First, we introduce the 'partial-wave' zeta functions for the modes inside the cylinder 
(which we shall call 'int' region) 



C^) = Ei m ;. ff^OO = £ *£» > far Re a >1, (13) 



n=l rt=l 



where j m ,n is the nth nonvanishing zero of J m and k m ^ n means the nth solution of eq. ([12] ), 
for a given m. The superscripts T> and M refer to Dirichlet and Neumann b.c. The analogous 
functions for the external modes ('ext' region) will be called (^(s), Cf'^O)- Looking at the 
whole problem, we have to take into account the angular degeneracies 

, 1, for m — 0, 
d(m) = { (14) 
2, for m > 0. 
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With them, we define the complete zeta functions for the circle: 

oo oo oo 



-mtJft 



m=0 n=l m=0 

oo oo oo 

m=0 n=l m=0 



(15) 



and 



oo 



(16) 



m=0 

With these ingredients, one builds up the total complete zeta function for the circle, i.e. the 
one entering eq.(^[), which is 



By the methods of ref.|| the functions (|I~3D have the following integral representations 
valid for — 1 < Re s < 0: 

a jr e /"OO r 

C art ' c ° nd (s) = - sin — / dx x- 8 - 1 In L part ' cond (m, a 
n 2 Jo L 



1* 



where 



part G { int, ext }, cond G {P,AT}-, 

L int ^(m,x) = v / 27rxe- a: J m (x), L mi ^{m, x) = V2^e' x I' m (x), 

/27 /27 ( 19 ) 

L c ^(m,x) = J—e x K m (x), L ext ^(m, x) = -J^ e x K' m (x). 

V 7T V 7r 

This is not enough yet, as we must go further to the left, until s = — 2. In order to extend 
the domain of validity of these formulas, it suffices to perform an adequate subtraction in the 
integrand, with the aim of separating the part that causes the divergent behaviour which stops 
us from reaching s = —2. One of the simplest possibilities is to add and subtract functions of 
the type 

x-^iy^m^x) + y 2 (m)t 2 (x)}, t{x) = (1 + x 2 )' 1 ' 2 (20) 
Specifically, we take 

y 1 (m) = -- (to - -j , y 2 (m) = -- I to - - 

y x (to) = ~2 V 4/ ' 2 ' ^ = ill 

yt^ conA {m) = -X nt,cond (m), K Xt * C ° nd (^) = ^2 nt ' C ° nd M, cond G {D,A/"}, 

(21) 



6 



and write 



x 



dxx 



In 



/■part, cond 

LPart, cond (m,a;; 
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s) = — sin — 

7T 2 



■ypart, cond (m)t(x) _ -ypart, cond ^2^) 
+jr t, cond (m) l_ B _0 + ypart, cond (m) 1^ _^ | > 

part G { int, ext }, cond G {"Z^A/"}. 



(22) 



Given that all the represented integrands have now the asymptotic behaviour x s 1 • O (t 3 (x)), 
the integrals are finite at s = —2. At this point, their specific values do not matter, as they are 
eventually multiplied by the zero coming from sin Further, since lim^ sin 
—it, all the expressions have a finite limit at s = —2. For instance, 
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-- I TO" 
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/■int, A/7 



-2) = i (^ 2 



(23) 



The first confirms a previous result in ref.JT^], while the other is a new result. 



Dirichlet modes. The integral in eq. (^) converges for s = —2 and can be obtained 
numerically. In order to compute the complete ( functions, ( |T5D and (|i~6|), it proves convenient 
to perform further manipulations. This job has already been done in ref.0. Although the 
relevant ( functions of that paper were obtained starting from expresions that are valid for 
— 1 < Re s < 0, we can make use of their final expresion (eq.(3.22) of ref.|§ with a[ ,T> = 
—o 1 ^ = — 1 for the interior modes, and a[ I,v = o 1 ^^ = — 1 for the exterior modes) provided 
N > 2. Adding the contribution of the internal and external modes, we obtain 
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(24) 



(25) 
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In ( |2"4| ) we have chosen N — 4. For the coefficients Sf(s,m) we find 

«Sf (-2, 1) = 1.34757127 x 1(T 3 , «Sf(-2, 2) = 2.83382470 x 1(T 5 , Sf (-2, 3) = 2.6944637 x 1(T 6 , 
Sf (-2,4) = 4.9500306 x 10 -7 , Sf (-2,5) = 1.31795554 x 10 -7 , Sf (-2,6) = 4.4524310 x 10 -8 , 
Sf(-2,7) = 1.7751724 x 10" 8 , Sf(-2, 8) = 7.99494 x 10 -9 , ••• (27) 

When m is large enough, Sf(— 2,m) decreases as m -6 , and can be approximated by 

v 19 7649 192349 15293983 

4 (- > m ) ~ 8960m 6 ~ 4730880m 8 + 82001920m 10 ~~ 2788065280m 12 ' ^ ^ 

In order to construct the complete zeta function, we need the series 
^ x s 2 (6 + 5sV , x s 2 (2 + s)(176- 268s -452s 2 - 113s 3 ) , 

+- sin — £ m-^f(s,m), (29) 

* 2 m=l 

where the appearance of the Riemann zeta functions Cr comes from identifying sums of the type 

oo 

m~ z = (r(z). Observe that, strictly speaking, such identification is only valid if Re (z) > 1. 

m=l 

Given that we are analytically continuing expressions involving s, we have to suppose that we 
may temporarily take s large enough so that this process be valid, before letting s — > — 2 at 
the end. It appears also the series 

oo 

^ rn 2 Sf (-2, m) = 0.001500509798 (30) 

m=l 

Expression ( p4j ) is not valid for m = 0, since it was obtained from a rescaling x — > mx and 
application of uniform asymptotic expansions in mx. However, the m = contribution can be 
easily found by applying formulas (|22|) and preceding, which yield 

C?00 = Co nt ' V (s) + C^is) = ~ + - sin ^ r dxx- 1 -* { \n{2xI (x)K (x)} - ~ ' 

16 7T 2 Jo I 8 

Numerically, we find 

f°° dxx( ln{2xI (x)K (x)} - - — ^— \ = 0.01096298110873 (32) 
Jo I 8 1 + x 2 ) 

Putting everything in (^9|)-(^2p together, we obtain, for s ~ —2, 

oo 

a s Cc 00 = E = 0.007725967(s + 2) + ((s + 2) 2 ) , (33) 

m=0 



1 + X 2 

(31) 



and, therefore, 



0.000614794033 



(34) 



where e v means the part of e coming from the contribution of (j^ only, and can be interpreted 
as the Casimir energy linear density of a scalar field under Dirichlet b.c. on the cylindrical 
surface. As expected, ( [54]) is free of divergences, which provides us with a nontrivial check of 
our result. The result in (|34j ) (and (|44]) and (^) below as well) is the lowest order result, which 
we expect to be modified by higher order (quantum, finite temperature...) effects. 



Neumann modes. The Neumann modes admit a very similar treatment. For the equivalent 
of (M) we find 
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f (s, m) 
(35) 



= S#(s,m). (36) 



(37) 



In order to obtain the complete Neumann ( function we need the sum 
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As we content ourselves with N = 4, we will only need S^(— 2, m), After some calculations, we 
have obtained: 



<(-2,4) = 
'(-2,7) = 



-2.3057255 x 10~ 3 , 
-8.2371882 x 10~ 7 , 
-2.9468381 x 10~ 8 , 



-4.755908 x 10~ 5 , <S^(-2,3) 
-2, 5) = -2.1903952 x 10" 7 , Sf(-2, 6) 
Sf (-2,8) = -1.3267927 x 10~ 8 , 



(-2,2) = -4.755908 x 10~ 5 



-4.495047 x 10~ 6 , 
= -7.3944610 x 10~ 8 , 
(39) 



9 



For sufficiently large m, (— 2,m) can be approximated by 

> r , , 9 1657 494891 10472667 , , 

S (—2 m) = I I f40) 

4 V ' ; 2560m 6 675840m 8 147603456m 10 1394032640m 12 ' V ; 



J2 m 2 Sf(-2, m) = -0.00256191227 (41) 



which implies that 

m=l 

Again, d35|) is not valid for m = 0, which has to be obtained separately: 

3s s . /7ts\ r°° i „ f, r „ s „, , s , 3 1 



tf{s) = ^ + £ sin (^) jT x" 1 -^ { In [-2^(^(3)] + g^T^}' ( 42 ) 



which for s ~ —2 can be approximated by 

tf( s ) = — + - sin — [-0.475214928727027 + 0(s + 2)1 

3^3 \ 

= — + 0.475214928727027 J (a + 2) + O ( (s + 2) 2 

8 V16 / v 

Gathering everything together, we find that the contribution of the Neumann modes is 



(43) 



0.01417613719 



e" = 5 . (44) 

cr 

which, as expected, is free of divergences. Here e 1 ^ means the part of e coming from the 
contribution of (j^ alone. It can be also viewed as the Casimir energy linear density of a scalar 
field subject Neumann b.c. on the cylindrical surface. 

Result and discussion. The total electromagnetic energy per unit of length is given by the 
sum of ([34]) and (fill), which yields 



£ e.m. = t V + e AT = — — 

a 2 



0.013561343 



which confirms the result obtained in ref. W%, where a high-frequency cutoff procedure was 



used, as well as our new results ( [34] ) and (|44]). The result in ( ^5|) is negative and indicates an 
ensuing attractive force. However, only the Neumann modes give a contribution of the same 
sign as the net result. The part from the Dirichlet modes is two orders of magnitude smaller 
and has opposite sign. In view of the analogous problem for the circle in 2 + 1 dimensions [f|, 
where the e.m. modes reduce to a set of Neumann modes only, this dominance is perhaps not 
so surprising. Nevertheless, in that case the result was infinite before applying the principal 
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part prescription (by zeta-function regularization, the type of divergence involved is seen to 
disappear when the number of space dimensions is odd). 

Curiously enough, one of the initial motivations for calculating e e m ' was the conjecture 
f0| that, being an intermediate situation between the parallel plates (attractive) and the sphere 
(repulsive), the cylinder might turn out to have e e ' m ' = 0. What we see is that the two 
contributions to the e.m. energy do indeed have opposite signs, although they fail to cancel 
each other as the authors of that speculation might have hoped long ago. 

In ref.[|llj the scalar field under Neumann b.c. was separately studied for the spherical 
case. Using the figure in that paper (and the ones in previous works for the other cases with 
a sphere), we can compare the Casimir energy per area unit Ec/A for the sphere and for an 
infinite cylinder with the same radius (a): 



a 3 • E c /A 


"D-scalar 


Af- scalar 


e.m. 


cylinder 
sphere 


+0.000098 
+0.000224 


-0.002256 
-0.017808 


-0.002158 
+0.003689 



Observe that, for the sphere, the e.m. case is not T> + jV because, in the presence of that 
surface, it decomposes into Dirichlet plus some special form of Robin — not Neumann — modes. 
Actually, the rate between the D and M contributions to the cylinder is —0.043. . .. In the 
spherical case, the ratio of l T> to A/"' has the value of —0.012 . . ., which has the same sign and 
order of magnitude, but is nearly four times smaller. Note that, in the case of the parallel 
plates, this ratio is +1. Comparing the e.m. results alone, the sphere has almost twice the 
absolute value of the cylinder. At any rate, since the energy for the cylinder is lower, the 
Casimir effect might tend to deform conducting spherical bubbles into cylindrical tubes. In 
a 'foam-like' universe model, one could even imagine a bubble becoming very thin in some 
directions and very long in another, as was conjectured in ref. ||. 

Regarding the two parts as separate results for scalar fields under different b.c, we 
observe that a change in the nature of the conditions can produce a drastic alteration, in size 
and sign, of the vacuum energy. Thus, the Casimir effect for a scalar field under Dirichlet b.c. 
would tend to expand the cylinder, but, after changing them into Neumann b.c, the effect 
would tend to make the cylinder contract. This is, too, the tendency of the e.m. field, as was 



duly pointed out by the authors of ref.|T7|]- In the same paper, they also said: Still eluding us 
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is the physical intuition to predict Casimir attraction or repulsion. Seventeen years later, the 
truth of this remark seems to hold on. 
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